Global quantum correlations in finite-size spin chains 
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We perform an extensive study of the properties of global quantum correlations in finite-size one-dimensional 
quantum spin models at finite temperature. By adopting a recently proposed measure for global quantum corre- 
lations [C. C. Rulli, and M. S. Sarandy, Phys. Rev. A 84, 042109 (2011)], called global discord, we show that 
critical points can be neatly detected even for many-body systems that are not in their ground state. We con- 
sider the transverse Ising model, the cluster-Ising model where three-body couplings compete with an Ising-like 
interaction, and the nearest-neighbor XX Hamiltonian in transverse magnetic field. These models embody our 
canonical examples showing the sensitivity of global quantum discord close to criticality. For the Ising model, 
we find a universal scaling of global discord with the critical exponents pertaining to the Ising universality class. 
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PACS numbers: 

Entanglement and criticality in quantum many-body sys- 
tems have been shown to be strongly and intimately connected 
concepts |[EI2]- The body of work performed with the aim of 
grasping the implications that critical changes in the ground 
state of a given Hamiltonian model have for the sharing of en- 
tanglement by the parties of a quantum many -body systems is 
now quite substantial [3]. This has resulted in important pro- 
gresses made in our understanding of the interplay between 
critical phenomena of interacting many-body systems and the 
setting up of genuinely quantum features. In turn, such suc- 
cess has proven the effectiveness of the cross-fertilization of 
quantum statistical mechanics by techniques and interpreta- 
tions that are typical of quantum information theory. 

Yet, it has recently emerged that the way correlations of 
non-classical nature manifest themselves is not necessarily 
coincident with entanglement, and a much broader defini- 
tion of quantum correlations should be given fllHl- This is 
encompassed very effectively in the formulation of so-called 
quantum discord as a measure striving to capture the above- 
mentioned broadness of quantum correlations @. In analogy 
with the case of entanglement, the relation between quantum 
discord and the features of quantum many -body models is fun- 
damentally interesting for the understanding of the role that 
the settlement of quantumness of correlations play in deter- 
mining the critical properties of such models. A systematic 
analysis in this sense, which has only recently been consid- 
ered [0-Qj]], is thus highly desirable. This is even more impor- 
tant given that some of the investigations performed so far hint 
at the fact that quantum discord is more sensible than entan- 
glement in revealing quantum critical points, even for systems 
that are not at zero temperature lfl2ll . 

This is a particularly relevant result, whose validity should 
also be checked for models that are both finite sized and at fi- 
nite temperature. The motivations for such an endeavor stem 
from the fact that, likely, the properties of quantum many- 
body systems will be addressed experimentally in systems 
consisting of, for instance, cold atoms loaded in optical po- 



tentials or trapped ions, as in very recent ground-breaking ex- 
periments II3I - I17I1 . At variance, the studies performed so far 
have mostly dealt with systems at the thermodynamic limit. 
Moreover, intuitively, one would expect global measures of 
general quantum correlations to be well suited to reveal the 
subtle features at hand here, given that some of the critical 
changes occurring in the lowest-energy state of many-body 
systems truly involve (quasi-)long-range influences among the 
parties. Such an analysis is made very difficult, both at a theo- 
retical and computational levels, by the lack of unambiguous 
measures of multipartite entanglement in mixed states. 

In this paper, we study the relation between criticality and 
global quantum correlations in finite-size systems at non-zero 
temperature by using a measure of g lobal quantum correla- 
tions recently put forward in Ref. I118TI and employed by some 
of us in |@1 for a quantum many-body system at zero tem- 
perature. As canonical examples, we study one-dimensional 
models that are of genuine physical interest due to the non- 
trivial features of their phase diagrams, such as the transverse- 
field Ising model, the open-boundary XX model in transverse 
magnetic field 0911. and the so-called cluster-Ising model in- 
troduced in Ref. 1201] . The latter interpolates between the stan- 
dard antiferromagnetic Ising Hamiltonian and a topologically 
ordered cluster phase. Our study shows the ability of global 
discord to detect critical points. Moreover, for specific cases 
among the examples addressed in our work, we bring evidence 
of a finite-size scaling for global discord and its derivative 
that are closely related to the behavior of macroscopic fea- 
tures such as the magnetization. 

The rest of this manuscript is organized as follows. In 
Sec. UJ we start our study by introducing both quantum dis- 
cord and its global version. In Sec. [II] we then move to the 
description of a set of physically relevant interacting quantum 
many-body models that will be studied against the content of 
global quantum correlations of equilibrium states at tempera- 
ture T + and present our results. Finally, in Sec.|IIl]we draw 
up our conclusions and discuss a few open questions that are 
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left to a future addressing. 



I. TOOLS FOR QUANTIFYING QUANTUM 
CORRELATIONS 

In this Section we introduce the fundamental mathemati- 
cal tools used in our study. We recall the definition of global 
discord given in Ref. ||18l| . and present a more agile expres- 
sion for the case of multipartite qubit systems. For the sake 
of completeness we briefly review the original formulation of 
quantum discord valid in the bipartite scenario. 



A. Quantum Discord 

We begin by reminding that, as originally proposed in 01, 
quantum discord is linked to the discrepancy between two 
quantum extensions of the concept of conditional entropy that 
are classically equivalent [4]. Let us consider a bipartite sys- 
tem described by the density operator p A B with p A (pb) denot- 
ing the reduced state of system A (B). The total correlations 
between A and B are quantified by the mutual information 

I(Pab) = S(p A ) -S(p A \p B ), (1) 

where S(pa) = _ Tr[p A log 2 Pa] is the von Neumann entropy 
and S(pa\pb) = S(pab) - S(pb) is the conditional entropy. 
By using a measurement-based approach, a second defini- 
tion of conditional entropy can be formulated. The appli- 
cation of a local projective measurement, described by the 
set of projectors {ft^} on part B of the system, results in 
the conditional post-measurement density operator Pab\j = 

(l A ®fLpPAB(U®fl B )/Pj, where Pj = Tr[(l A ®ft^)p AB ] is the 
probability associated with the measurement outcome j. We 
can thus define the measurement-based conditional entropy 
S(pAB\n J B ) = XjPjS(pA\j) withp AL/ = Ti[n J B p AB ]/p j , which 
leads us to the so-called one-way classical information ||5| 

J(p A B) = S(pA)-S(pAB\tli). (2) 

The difference between quantum mutual information and clas- 
sical correlations, minimized over the whole set of orthogonal 
projective measurements performed on B, defines quantum 
discord as 

£> B ^ A (Pab) = inf IKpab) - J(Pab)1 (3) 

By noticing that the original definition of discord 0] 
can be rewritten in terms of relative entropy SifiiWpz) = 
Tr[pi lo g 7 p il - Tr[pi log 2 p2] between two generic states p\ 
and p2 flal and by symmetrizing its definition through the 
introduction of bilateral measurements ft^ ® ft* |2lll . we in- 
troduce 

^{PAB) = min [S(p AB \\Tl(p AB )] - V S(pj\\hj(pj)) (4) 
(n>n») hr B 



with tl{pAB) = Zj,k(^i ® % )p AB {tl i A ® ft*) Eq. © expresses 
discord as the difference between the content of quantum cor- 
relations ascribed to a multi-local measurement process and 
the sum of the relative entropies for each reduced state of the 
system. The minimisation is required, clearly, to remove any 
dependence on the local measurement bases. The absence of 
global quantum correlations would make Eq. <j4j identically 
null. 



B. Global Quantum Discord 

Eq. (0| is the starting point for the formulation of global 
discord (GD) d 

QDipj) = min js (prl|ft(pr)) - g S (p;||ft>;))J , (5) 

which quantifies the global content of non-classical corre- 
lations in the state pj of an A^-party system. Here pj = 
Tr' [pr] is the reduced state of qubit j [we use Tr' for the 
trace over all the qubits but the j ], fljipj) = 2/ft'pjft'., 

ft(pr) = Zk tl k p T tl k , tl k = S^ftf , and k stands for the string 
of indices {k\ . . . k^). The minimization inherent in Eq. (0 
is performed over all possible multi-local projectors ft*. In 
Ref. lfl8ll it is shown that QD(pr) > 0, its maximum value 
depending on the dimension of the total Hilbert space at hand. 
Recently, a monogamy relation relating global quantum dis- 
cord in a multipartite setting and pairwise correlations evalu- 
ated by quantum discord has been introduced in Ref. 1E2I . 

The explicit computation of the formula in Eq. (0 is in 
general a difficult problem. However, the task can be greatly 
simplified by writing the multi-qubit projective operators as 
ft* = 7? |k) (k| 7?' . Here (|k)} are separable eigenstates of 
= ®JLjcH with & q . the q = x,y, z Pauli operator, and 

'R is a local multi-qubit rotation % = ® N . =l Rj{6j,(f)j) with 

Rj(9j,(f>j) = cosOji + i sin 6j cos <f>j& y + i sin 9j sin <pjfr x the 
rotation operator (of angles 9 and (f>j) acting on the y'-th qubit. 
Analogously, the set of local projective operators on the y-th 
qubit is written as ft^ = Rj \l) (l\ £t with \l) (1-0,1) the eigen- 
states of and where, for convenience, we have dropped the 
dependence of the rotation operators on their respective an- 
gles. As shown in some details in the Appendix, the introduc- 
tion of these quantities allows us to reformulate GD as 

I N 1 2^-1 "I 

§D{p T ) = min i P" l«g 2 P" ~ 2 Pt l°g 2 Pr 

m U =I '=° *=° i (6) 

N 

+ J]s(pj)-S(p T ) 

7=1 

withpf = (kl'RVr^lk) andp^' = </|^jp^|Z). Despite its 
involved form, Eq. (|6]l greatly reduces the computational ef- 
forts needed to evaluate QT)(pt)- 
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II. QUANTUM CORRELATIONS AND CRITICALITY IN 
QUANTUM SPIN CHAINS 

To examine the qualitative and quantitative features of 
quantum correlations in spin systems, we will focus on three 
different models: the Ising, cluster-Ising and the XX models. 
Here we are interested not only in the zero temperature case, 
but also in exploring the thermal effects on such finite sized 
quantum systems that exhibit critical behavior. For the size of 
the systems we will be considering, the explicit thermal state 
can be directly calculated via its canonical ensemble and is 
given by the Gibbs state [throughout this manuscript we take 
units such that % = ks = 1]: 

-«/r 

Q(T) = — (7) 

with ( h the Hamiltonian describing the interaction, T the ef- 
fective temperature, and Z = Ti[e^^ T ] the partition function. 

A. Transverse field Ising Model 
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FIG. 1 : QT) for the Ising model at zero temperature. From bottom 
to top curve, L goes from 3 to 1 1 spins. At B = the ground state 
of the spin model is a GHZ (corresponding to the ferromagnetic case 
with no symmetry breaking), thus giving QT) = 1 lfl8h . In the para- 
magnetic configuration (B » J), QT) goes to zero together with any 
non trivial spin correlation. 



We start our analysis considering the quantum Ising model 
in the zero-temperature case. The behavior of bipartite and 
global correlations in the transverse spin-1 /2 Ising model has 
attracted considerable interest so far. Entanglement (jH l23ll . 
non-locality {24}], and bipartite quantum discord have been 
studied for this model. More recently, the scaling of entangle- 
ment spectrum of a finite-size spin- 1/2 Ising chain near its 
critical point has been studied 112511 . 

Here, in line with some of the studies mentioned above, we 
shall consider a one dimensional system with periodic bound- 
ary conditions. The Hamiltonian for a chain of L spin- 1/2 
particles reads 

L L 

n--jYj & ^ + B Tj^ (8) 

1=1 !=1 

with the condition L + 1 = 1. In the limit B/J — > 0, the 
ground state of this model is locally equivalent to an L-spin 
GHZ state I23I1 . As B increases, the entanglement in the 
ground state soon disappears, as the spins tend to align along 
the direction set by the magnetic field and, in the thermody- 
namic limit, the system undergoes a quantum phase transition 
at B/J = 1. The nonlocal nature of the quantum correlations 
within the ground state of this model has been studied in ll24Tl 
and found to be extremely sensitive to temperature, a feature 
shared with the entanglement. Needless to say, this does not 
imply that all non-classical features in the correlations shared 
among all the spins disappear with temperature and we shall 
demonstrate that global quantum discord is indeed able to sig- 
nal the structural changes in the sharing of quantum correla- 
tions even at T + 0. 

Technically, the evaluation of GD for the Ising model en- 
joying the symmetries mentioned above offers room for a few 
interesting considerations. First, as the model in Eq. (|8) is 



characterized by a real symmetric Hamiltonian, GD is com- 
pletely independent on the set of <pj angles, which do not play 
any role in the minimization necessary to calculate Eq. ((5). 
Second, as we have taken periodic boundary conditions, the 
system is translational invariant. Such invariance has conse- 
quences on the relation among Qj, making GD invariant under 
cyclic permutations of such variables. Finally, we have gath- 
ered numerical evidence of even a higher degree of symmetry 
of the GD function, at all lengths and values of the temper- 
ature, in that the optimal dj all take the same value 9, which 
depends on the magnetic field. We have 6 = [n/4] for small 
[large] values of B. The transition between these two values 
is sharp and happens in proximity of the critical point, thus 
showing how the changes induced at criticality are reflected 
in the structure of GD. We have run our simulations by min- 
imizing over all the possible different angle configurations, 
finding perfect agreement with the results corresponding to 
this explicit choice of angles. 

Fig. Q] shows the amount of quantum correlations, quanti- 
fied by GD, as a function of the ratio B/J between the mag- 
netic field intensity and the Ising interaction constant. We 
study rings with L = 3, . . . , 1 1 whose GD curves share the 
same value at zero magnetic field. This agrees with the ground 
state being an L-spin GHZ, for which QT) = 1 regardless of 
the system size. As B/J tends to 1 the global discord increases 
reaching a maximum at different positions depending on the 
length of the chain. In the paramagnetic phase achieved for 
B » /, all the spins align along the direction of the magnetic 
field, so that all quantum correlations disappear. At smaller 
values of B, however, the global sharing of quantum corre- 
lations undergoes significant changes which result in the ap- 
pearance of a maximum, whose height and position is a clear 
function of the size of the system. 

The core part of our analysis consists of the study of the 
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changes in the behaviour of GD for rings prepared in thermal 
states. We consider two different cases with effective temper- 
ature equal to T = 0.05 [cf. Fig. 12(a)] and T = 0.1 [cf. 
Fig. |2](b)]. At non-zero temperature, the quantum correla- 
tions that are present in the ground state at B = are de- 
stroyed and QT)(pt) = VL. Overall the height of the curves 
decreases with increasing temperature. At low temperatures 
[cf. Fig. [2] (a)], the position of the maxima of QT>(pj) is ex- 
tremely close to those at T = 0, while higher temperatures 
induce a shift in the maxima of each curve [the effect is al- 
ready visible in Fig.[2](b)]. The optimal angles for the GD of 
thermal states are the same as those for T — 0. 

In line with the studies performed on the scaling of entan- 
glement in quantum spin systems lErjll . it is interesting to study 
the way global quantum correlations scale against the number 
of elements of the multipartite systems that we are addressing 
here. To this aim, we perform a finite-size scaling analysis of 
GD in the proximity of the critical point B,„/J for a finite-size 
transverse Ising model at T — 0. We thus study the derivative 
of Q"D with respect to B finding that, in proximity of the finite- 
size critical point, it is a function of U(B - B„,)/J satisfying 
the scaling ansatz ll8l l27il : 
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L-"f[U(B-B,„)/J] 
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FIG. 2: QD for Ising rings containing from 3 to 10 spins at non-zero 
temperature: (a) T = 0.05, (b) T = 0. 1. GD is null at zero magnetic 
field. For increasing temperatures the maximum values of the curves 
decrease. 
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FIG. 3: Derivative of GD with respect to the magnetic field for the 
Ising model at zero temperature. The data points (for L = 4, . . . , 11) 
are scaled according to the number of spins of the rings. B m is the 
critical value of the magnetic field for a finite-size chain of length L. 
Close to the critical point the curves at all values of L collapse to a 
common function, witnessing universality. 



where v — 1 is the correlation length divergence critical expo- 
nent of the corresponding Ising universality class. We are able 
to make data collapse already for very small lengths (L < 11) 
and we obtain w*-1.5 while f(x) is approximately quadratic 
for x ~ [cf. Fig. [3). This is quite a remarkable result, as it 
shows that GD scales with universal critical exponents in the 
proximity of a quantum phase transition. 



B. Cluster-Ising model 

We now consider a model, recently proposed and studied 
in Ref. lEoTl . which combines competing effects coming from 
an antiferromagnetic Ising and a three-body cluster interaction 
according to the Hamiltonian (with periodic boundary condi- 
tions) 



CI 



';+] 



+ A 



(10) 



The three-body term in Eq. (fTOl is responsible for the set- 
ting of long-range entanglement, which has been recently 
related to topologically ordered states 112 811 . while the sec- 
ond term tends to localise entanglement to nearest-neighbor 
pairs of spins. Such competition makes the model undergo a 
second-order quantum phase transition at J/ A = 1 with the 
ground state of the system passing from an Ising antiferro- 
magnetic phase (for J/ A «: 1) to a cluster-like one (achieved 
at J/ A » 1) endowed, as said above, with long-range en- 
tanglement and topological order. Such transition, which is 
not in the Ising universality class, has been characterized in 
Ref. I20I l29ll by means of a global geometric measure of en- 
tanglement. 

An interesting point to notice is that, for a cluster-Ising 
model at non-zero temperature, neither the two-spin nor the 
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multipartite entanglement (as measured by tangle) is able to 
signal the quantum phase transition, as they are either identi- 
cally null (two-spin entanglement) or equal to a constant (mul- 
tipartite tangle). Here we will make use of global quantum 
discord to study the occurrence of critical structural changes 
in the correlation-sharing structure of the model even at fi- 
nite temperature, showing the effectiveness of GD in the task 
of revealing such modifications at criticality. Technically, the 
problem of finding the minimum in Eq. (0 is more difficult 
to tackle than in the Ising case. In fact, the model is charac- 
terized by a lower degree of symmetry (due to the presence 
of all the Pauli spin operators in Eq. (fTUl i) which forces us to 
minimize Eq. © using all angles {#,} and \<pj). For J/ A <K 1 
(Ising phase) the optimal angles are {6j = n/4} and {<pj = n/2) 
analogously to the ferromagnetic phase of the Ising model (the 
value of the phases <pj is due to the different Ising coupling of 
this model). In the opposite regime J/ A » 1 (cluster phase) 
the optimising angles depend on the number of spins. 

We start the description of our result by analyzing Fig. |4] 
where we plot the GD as a function of the ratio J/ A for 
L = 3, . . . , 10 and at T — 0. In the limit of vanishing cluster- 
like contribution, the model contains the Ising two-body inter- 
action and we correspondingly recover the results obtained in 
Sec. Ill Al for zero transverse magnetic field: the ground state 
is locally equivalent to an L-spin GHZ state and QT) = 1 re- 
gardless of the size of the chain. In the opposite asymptotic 
regime, where the three-body model dominates over the Ising 
term, the ground states are cluster states on an L-site ring lat- 
tice [f3Qh - Even at moderately large choices of J /A, the value 
taken by GD agrees very well with the expectations for size- 
L cluster states, which are equal to 1 , 2, 3, 3, 4, 4 for L going 
from 3 to 8, respectively. Away from such limits, GD behaves 
in a peculiar way with the size of the system. Most of the 
cases that we have considered in our analysis display a non- 
monotonic behavior with a peak occurring in proximity of the 
critical point. However, among the values of L considered 
in our calculations, the cases of L = 3 and L — 6,9 behave 
differently, with the GD being practically constant or with a 
sharper maximum at values of J/ A significantly away from 
the critical point. Moreover, these special cases give rise to a 
few crossings with the curves associated with both lower and 
larger rings (for instance, the curve corresponding to L = 6 
crosses both those for 5 spins and 7 spins). 

We believe that the occurrence of such pathological behav- 
ior when L is a multiple of 3 should not be regarded as ac- 
cidental, but rather as a signature of the distinctive features 
of the cluster-Ising model for these lengths of the system. In 
fact, the study (at finite size) conducted in Ref. [29] shows 
that, differently from the thermodynamic limit, the x- and z- 
correlations in the model, as well as the magnetization along 
Z axis, vanish for L that is a multiple of 3 (Ref. l29ll discusses 
explicitly the case of L = 6 and 12). Although, given the com- 
plexity of the task, it is implausible to formulate an analytic 
expression of GD from which the role played by such corre- 
lations can be clearly extracted, we conjecture their relevance 
in the determination of the functional form of GD at moder- 
ate values of J /A, where the differences with respect to any 
other system-size appear to be more striking. In turn, given 
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FIG. 4: QT) for the cluster-Ising model at zero temperature. We took 
L = 3, . . . , 10, growing from the bottom to the top curve as they 
appear in the J/ A » 1 part of the plot. In the second case QT) for 
cluster states with an increasing number of spins starting from 3 is 
equal to 1, 2, 3, 3, 4, 4. 



the agreement between the calculated GD and the predictions 
valid for the ground state of the Ising and cluster models at all 
values of L, this analysis suggests that the above mentioned 
correlators are not heavily relevant for the calculation of GD 
deep in both such phases. While this observation could help 
in studying this figure of merit analytically, a less phenomeno- 
logical approach to this interesting points goes beyond the 
scopes of this work and remains to be addressed in further 
studies on this matter. 

We conclude our analysis of this model by addressing now 
the case of T + 0, as done in Fig. [5] While GD vanishes in the 
antiferromagnetic phase, it persists to the effects of temper- 
ature in the cluster one. This is in line with what has been 
found for thermal cluster states, for which an L-dependent 
critical temperature exists, below [above] which distillable 
and long-range [bound] entanglement is found in the cluster 
state lHHm. 

Even more strikingly, though, the structure observed at 
T — survives, qualitatively unaltered, at T + 0. Actually, the 
quantitative differences between the results associated with 
the L-spin ground state and the corresponding thermal equi- 
librium state are negligible [in terms of both the position of 
the maximum of GD on the J/ A axis and their actual value, 
cf. Fig. [5] (a) and (b)] even at values of T for which the GD 
of an Ising chain was found, in Sec. Ill Al to be sensibly dif- 
ferent from that of the T = case. Needless to say, at much 
larger values of the temperature the peaks close to the critical 
point are smeared out into a broad and flat GD curve. This 
demonstrates the claimed effectiveness of the figure of merit 
addressed here in signalling the effects of criticality on the 
sharing of quantum correlations in thermal-equilibrium states, 
thus reinforcing the significance of the analysis conducted so 
far along the lines of combining quantum many-body physics 
and discord-related quantifiers []7|-|9|]. 
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C. Open-chain XX model 

Finally, we address the case of an open-ended chain inter- 
acting via an XX term in the presence of a transverse magnetic 
field 



n 



XX 



-ilf 



(T; cr. 



+ cr. cr, 



(11) 



As studied in Ref. 113311 and shown for L — 8 in Fig. [6] the 
energy spectrum of this model is quite rich and encompasses 
quite an interesting series of crossings among its energy eigen- 
states occurring as the ratio B/J is varied. Correspondingly, 
the ground state of the system changes and can be classified 
in terms of the number of such level crossings. In the thermo- 
dynamic limit L — > oo, 'Hxx exhibits a Berezinskii-Kosterlitz- 
Thouless (BKT) transition at B/J = 1 from a critical phase, 
characterized by in-plane quasi-long range order, to a para- 
magnetic one. Correspondingly, two-spin quantum entangle- 
ment (as measured by concurrence) reduces as the distance be- 
tween the spins is increased, and vanishes at the critical point, 
leaving a fully factorized ground state. 

Fig. [7] shows GD as a function of the external magnetic 
field over the two-body interaction constant. As discussed 
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FIG. 5: QT) for the cluster-Ising model at non-zero temperature. We 
took L = 3, . . . , 10, growing from the bottom to the top curve as they 
appear in the J/A » 1 part of the plot. In panel (a) [(b)] we took 
7 = 0.05 [J = 0.1]. 
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FIG. 6: Spectrum of the model defined in Eq. i ll It against the ratio 
B/J for L = 8 spins. Dashed vertical lines identifies the points of 
occurrence of level crossings in the model. 



above, for high magnetic field the system is in the paramag- 
netic phase and no correlations are present. For low magnetic 
field the GD displays a step-wise behavior, jumps occurring in 
correspondence of the level-crossings that redefine the ground 
state of the system (which are evident from the spectrum of 
the model). That is, GD tracks the structural changes in the 
ground state of the spin system as B/J varies. Moreover, as L 
grows, GD goes to zero at values of B/J increasingly closer 
to the critical point. Already for the modest value of L = 9, 
the difference between the value of B/J at which QT) = and 
B/J — 1 is only 5%. A non-zero temperature smoothens the 
sharpness of the jumps occurring in GD and reduces its am- 
plitude as shown in Fig. [8] Yet, the series of level crossings 
at which the ground state changes, as well as the BKT point 
are still clearly visible as dips between quasi-plateaux and a 
gradual yet quick decrease of QT), respectively, thus proving 
the effectiveness of GD as a figure of merit for signalling crit- 
icality at T + 0. It should be noted that, for B = T = 0, the 
ground state is doubly degenerate. In order to evaluate GD 
we have thus taken the linear combination of such degenerate 
states that smoothly provides the values associated to T + 0, 
where the degeneracy is lifted. 
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FIG. 7: QT) for the XX model at T = 0. We have taken L = 3, .... 9 
(from bottom to top curve). 
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III. CONCLUSIONS 

We investigated the behavior of a measure of global quan- 
tum correlations in three finite-size one-dimensional quantum 
many-body systems close to their critical points. We have ma- 
nipulated the expression for global quantum discord so as to 
adapt it to the case of quantum spin chains. This has allowed 
us to study the global quantum discord in thermal states of 
moderately large quantum spin models, demonstrating its ef- 
fectiveness in spotting the critical changes in the state of the 
system as a function of relevant parameters. Furthermore, for 
the Ising model we have been able to put forward evidences of 
a finite-size scaling behavior characterized by universal criti- 
cal exponents of the Ising universality class. For the cluster- 
Ising model, GD provides an alternative and powerful tool for 
the signalling of criticality superior to entanglement measures 
in both reduced and global forms. Finally, for an open XX 
chain, we have been able to track the discrete number of struc- 
tural changes for the ground state as the transverse magnetic 
field varies. 

Our analysis sheds light also on technical aspects related 
to the calculation of global discord in multipartite spin sys- 



tems enjoying some degree of symmetry. For the transverse 
Ising model with periodic boundary conditions, we have pro- 
vided strong numerical evidence that identical local projec- 
tions should be implemented in order to attain the global max- 
imum inherent in the definition of GD. Moreover, the az- 
imuthal angles </>j are shown to be immaterial for this task. 
Differently, for the cluster-Ising model in Eq. ([TOt , deep in the 
cluster phase, more complicated combinations of minimizing 
angles are found. Although, it would be interesting to find a 
relation between the symmetries of the model in considera- 
tion and the angles minimizing the global discord, a compre- 
hensive solution seems highly non-trivial and goes beyond the 
scope of this study. 
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Appendix 

Here we present the main steps of the derivation of the ex- 
pression of the global discord for N spin-system presented in 
Eq. ©. The starting point is writing the multi-local projectors 
as fl k = "R |k) (k| ft where we remind that j|k)| are multi-local 
(separable) eigenstates of the operator £- = ®^1 1 <7~ and "R is 

a multi-local rotation R = ®^ =1 Rj- With this definition the 
terms of the relative entropy of the form Tr[pi log 2 p2] (we 
call it here mixed terms) in Eq. ((5} can be rewritten as 



Tr [p T log 2 ft(p r )] = Tr 



Tr 



p r log 2 £ft|k><k|&Vr&|k><k|& 1 
p T log 2 R J] (lk>P?<kl)^ 

k 

Pr2 lo g2P"l k >< k l 



Tr 



( 12 > 

where we define p T = ftp T R, pf = <k| ft p T fo |k>, we use 
that if RAft = A then Rf (A) ft = f (A) and the cyclic prop- 
erty of the trace. Using exactly the same line of reasoning for 
the mixed terms of the relative entropy for the j-th qubit we 
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find that 
Tr[p;log 2 n;(p;)] = Tr 



Pj log 2 J]Rj\l)(l\R]pjRj\l)(l\R] 



/=() 



(13) 



where p u . = {l\ R'jPjRj \l) and |/) being the two eigenstates of 
cH. By putting all the terms together and taking into account 
the minimization, we obtain the expression for global discord 
given in Eq. ©. 
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